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We estimate final state interactions in the B-meson decays into two pions by the Regge model. We 
consider Pomeron exchange and the leading Regge trajectories that can relate intermediate particles 
to the final state. In some cases, most notably B and B tv'^tt , the effect is relevant and 

produces a better agreement between theory and experiment. 

[/^'. PACS numbers: 13.25.Hw 
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C\l . I. INTRODUCTION 

^ ■ Measuring the angle a of the unitarity triangle is one of the major challenges of the _B-factories BaBar at SLAG, 
BELLE at KEK and the future LHG at GERN. B — > tttt decay channels were identified long ago as a promising 
candidate for the extraction of the angle a. Though other channels were subsequently investigated, the B decay into 
two pions is still object of intense studies, both experimental and theoretical. The task of determining precisely the 
I ,. angle a is complicated by the problem of disentangling two different hadronic matrix elements, each one carrying its 
^ ■ own weak phase. They are usually referred to as the tree and penguin contributions. 

The theoretical uncertainty in the evaluation of these terms stems from the approximate schemes used to compute 
OO the relevant four-quark operators between the hadronic states. In principle one could avoid these uncertainties, 
_ fixing all the hadronic parameters by a simultaneous measurement of physical observables. This problem has been 
addressed by several authors; for example in Ref. [1] the use of isospin symmetry among the various B — > tttt channels 
was envisaged. This program represents a significant experimental challenge and is therefore useful to have some 
, theoretical indications on the results. In the Standard Model one expects the dominance of tree diagrams in B — > tttt 
' decays, differently from the B ttK decay channels, where penguin contributions should play a key-role. Thus one 
naively expects that for B tttt the hierarchical structure follows the analogous hierarchy of the Wilson coefficients, 
namely 

6(7r"^0)<6(7r+7r-)w2 6(^+7r") . (1) 
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As discussed in section IV below, present experimental data are at odds with Eq. (1). There can be several factors 
. !Z| ' leading to violations of the expectation (1). First of all the role of penguin operators should not be neglected. Second, 
■ one has to go beyond naive factorization and use more sophisticated schemes taking into account QGD in factorization, 
?H ' for example the BBNS approach [2, 3] or the Soft-GoUinear-Effective- Theory (SGET) [4-7] (for a recent discussion of 
B decay into two light mesons in the framework of the SGET see [8]). Finally, (1) does not take into accoimt final 
state interactions (FSI). To this issue the present paper is devoted. FSI are long distance effects that in some cases 
might play a significant role; for example in B decays into two light mesons a source of long-distance contributions 
is provided by the charming penguin diagrams that might produce the discrepancy between experimental data and 
the naive factorization findings. Gharming penguins are contributions where the final state is formed only as an 
effect of a rescattering process, and is preceded by the formation of an intermediate state containing a cc pair [9-14]. 
In particular for decay channels with a strange light meson in the final state, e.g. B ttK, these long-distance 
contributions are not numerically suppressed. In fact, the Gabibbo-Kobayashi-Maskawa (GKM) matrix elements 
produce an enhancement ~ jVcbT^sl/IKjbl^*,! which can compensate the parametric suppression predicted by QGD 
factorization. 

The role of charming penguins in _B — > tttt is less clear. Due to the lack of the above-mentioned enhancement, on 
general ground one expects a minor role in the B — > tttt decay modes. On the other hand in [8] their role is found to 
be significant. This matter should be settled, but in any case FSI must be taken into account, be they dominated by 
charming penguins or by other rescattering processes, involving non-charmed particles in the intermediate state. 

The most accurate way to take into account FSI in hadronic B decays is provided by the Regge model of high 
energy scattering processes, which can be applied to hadronic B decays due to the rather large value of s = m^. The 
advantage of the Regge approach is to evaluate the rescattering not by a Feynman diagram, but by unitarity diagrams 
and the Watson's theorem [15]. In particular there is no extrapolation of low energy effective theories to the hard 
momenta regime and therefore no need to introduce arbitrary cutoffs in the light meson momenta, because the high 
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energy behavior is completely under control. Some studies on the application of the Regge model to B decays are 
in Refs. [16-18]. Elastic contributions to high energy scattering are dominated by the Pomeron exchange, while the 
inelastic channels get contributions from both Pomeron and Regge trajectories. Also charming penguins find a place 
in this scheme, provided one introduces also charmed Regge trajectories, as for example in the study performed in 
[19] for the charmless B decay into two light vector mesons. 

The aim of this letter is to extend the results of the Regge model to the B ^ tttt decay modes. We will show that 
there are indeed significant rescattering effects in the B 7r°7r° channel, a decay mode that is suppressed in naive 
factorization. We include several intermediate states: tttt, pp, aiw, DD and we find that the largest contribution 
comes from the tttt and pp intermediate states with p and a2 Regge exchanges respectively. On the other hand we find 
no significant role of charming penguin contributions. The suppression of charming penguins in comparison to other 
terms is produced because Regge charmed trajectories have a negative intercept a(0) and therefore a suppression 
factor (s/so)"(°) (so ~ 1 GeV^, a threshold). 

The plan of the paper is as follows. In section II we evaluate in the factorization approximation the bare amplitudes, 
including tree and penguin contributions with no FSI (in particular no charming penguins). In section III we discuss 
rescattering effects parameterized by the Regge model. Finally, in section IV we present our numerical results and 
discuss them. 



II. BARE AMPLITUDES 

The TTTT final state can be reached from several intermediate states via rescattering. Clearly one should select the 
most prominent channels. Among the inelastic channels we single out the decays B ^ pp and B — > aiTr since they 
have large branching ratios. For example: B{B+ -> p"p+) = (26.4±6.1) x 10"*^; B{B" p- p+) = (30.0±6.0) x 10"^ 
B{B^ ^r^"^) = (42.6 ± 5.9) x 10~^ [20]. To these decay processes we have to add the elastic B ^ nn channels, 
though they have smaller branching ratios. We also add the D'^*^D^*\ having in mind a discussion on the charming 
penguins. In conclusion the final state interactions that we consider are the elastic scattering tttt — > tttt, and the 
pp TTTT and aiTT tttt and the £)(*)£)(*) tttt inelastic channels. 

We evaluate bare amplitudes in the factorization approximation. To do that one needs different input parameters. 
The non-leptonic hamiltonian is well known and we do not repeat it here, see e.g. [21]. For the Wilson coefficients we 
use: ai = 1.029, 03 = 0.140, and (03, 04, ag, ag, a-j, as, ag, aw) = (33.33, -246.66, -10, -300, 1.95, 4.81, -93.30,-12.63) 
xlO"^ [22]. We use for the parameterization of the CKM matrix sin 6*12 = 0.2243, sin 6*23 = 0.0413, sin0i3 = 0.0037 
and ^13 = 7 = 1.05 [23]. As for the form factors and constant decay we use /jr = 0.132 GeV, fp = 0.210 GeV, fai ~ 0.21 
GeV (see the discussion in [24]) F^^^iO) = 0.26, A^^''{0) = 0.26, a!^^'\0) = 0.23 , V^f ^"^(0) = ^^^''(0) = 0.39, 
where we use the notations of [25] for the B p transition and the parameterization of Ref. [26] for the B ain 
matrix element. All the other parameters are taken from [23]. We get in this way the results of Tables I-III (notice 
that units of Tables I-II are 10"^ GeV, those of III are 10"'^ GeV). 



Process 


At 


Process 


At (A = +l) 


Ab (A = -1) 


Ab (A = 0) 


B+ ^ tt+tt" 
B° ^ 77+7r- 


+2.02 - 1.24 i 
-0.41 + 0.053 i 
+2.43- 1.74 i 


B° ^ p°p° 
B°^p+p- 


-0.02 + 0.01 i 
+0.004 - 0.001 i 
-0.02 + 0.02 i 


-1.1 + 0.65i 
+0.20- 0.07 i 
-1.31 +0.85 i 


+4.49 - 2.76 i 
-0.83 + 0.31 i 
+5.53 -3.59 i 



TABLE I: Bare amplitudes for B ttiv and B — » pp. Results in 10 GeV; A = ±1, refers to the helicities of the vector 
particles. 



Process 


Ab 


B+ ^ a+7r" 


+3.4 - 2.0i 


B+ ^ aU+ 


+2.2 - 1.6i 


B° ^ a?7r° 


-0.60 + 0.20 i 


5° ^ a+n- 


+4.2-2.7i 


B° ^ a^7r+ 


+3.4- 2.4 i 



TABLE II: Bare amplitudes for B — > ai7r; ai with longitudinal polarization. Units are 10 ^ GeV. 
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Process 




Process 


At 


B+ ^ D+D" 


-2.8 1 


B+ 


^ D'+D" 


+2.5i 


B° D+D- 


-2.8i 




D*+D- 


+2.5i 


B+ ^D'+D'" 


-2.9i 


B+ 


^ D+D*'> 


+2.5 i 




-2.9 i 




D+D*- 


+2.5i 



TABLE III: Bare amplitudes B — > D^*'/)'*'; vector particles have longitudinal polarization. Results in 10 GeV. 



III. FINAL STATE INTERACTIONS AND REGGE BEHAVIOR 

Corrections to the bare amplitudes due to final state interactions are taken into account by means of the Watson's 
theorem [15]: 

A = VSAb (2) 

where S is the S'-matrix, Ab and A are the bare and the full amplitudes. An application of the Watson's theorem was 
first discussed in [17] and subsequently applied to other decay channels in [18] and [19]. We briefly review here the 
formalism. 

The two-body S'-matrix elements are given by 

sip =Sij + 2i^Tl'\s) , (3) 

where i,j run over all the channels involved in the final state interactions. The J = 0, isospin I amplitude Tlp{s) is 
obtained by projecting the J = angular momentum out of the amplitude T^p{s,t): 

t}P{s) = -^^ I dtT};>{s,t) . (4) 



['dtTl'\s,t) . 



pj ,£j and t± are defined in Rcf. [19]. For the channel B — > tttt we only have the / = and 7 = 2 transition 
amplitudes; the decay amplitude B — > tt+tt" is only 1 = 2. 

The phenomenological basis for the application of the Regge model of final state interactions is the large value of 
s = m^; therefore a Regge approximation based on Pomcron exchange and the first leading trajectories should be 
adequate. The Pomeron term contributes to the elastic channels. As discussed in previous section for the inelastic 
case we include only channels whose bare amplitudes are prominent. 

In conclusion in the present approximation wc will include, besides the Pomcron, the p and 02 (almost) exchange- 
degenerate trajectories and tt Regge trajectories. We shall discuss in subsection IIIB the role of charmed Regge 
trajectories in parameterizing charming penguins. 

For the Pomeron contribution we write (neglecting light meson masses) 

5=l + 2iT^(s), T^(s) = _i_y°T^(s,f)dt, (5) 

and we use the following parameterization [18, 27]: 



T^{s,t) = - (i-^g{t) e-'('^/2)«H*) , (0) 

with So = 1 GeV^ and ap{t) = 1.08 -I- 0.25t {t in GeV^), as given by fits to hadron-hadron scattering total cross 
sections. For the Pomeron residue f^^ we assume factorization with a t-dependence given by [18, 27] 

The additive quark counting rule allows to compute the Pomeron-pion residue in terms of the Pomeron-nucleon ones. 
This gives [16, 18): 

- - 5.1 . (8) 
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As observed in [17] inelasticity effects play an important role in the determination of the FSI phases. Parameterizing 
them as in Ref. [17] by one effective state, with no extra phases would allow to write the S'-matrix as follows (neglecting 
a small phase = —0.01 in \/l + 2iT'^): 

^ ^ ' 0.82 i 0.62 



0.79 0.64(1 + i) 
0.64(1 + i) 0.79 



(9) 



This shows that even neglecting the effect of the non leading Regge trajectories, final state interactions due to 
inelastic eflfects parameterized by the Pomeron exchange can produce sizeable strong phases. This result agrees with 
the analogous findings of Rcfs. [17] and [18]. However this method is not useful to evaluate rescattering effect in weak 
decays. Therefore we prefer to parameterize inelastic effects by Regge trajectories. 



A. Regge trajectories 

Let us now consider the contribution of the leading Regge trajectories. Including Regge trajectories the S matrix 
can be written for the generic B ^ nn case as follows 



(10) 



Here P indicates the Pomeron contribution discussed above. Since the Pomeron is much larger than the others we 
make the approximation 



A{B ^ 7r7r)(^) « Vl + 2iT'PAl^^ + 



1 



2^1 + 2iT'P 



(11) 



k n 



Here the sum over k refers to the various intermediate states contributing to the final state tttt; I is the isospin index. 
We write the Regge amplitudes as follows {TZ = p ,a2 
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We assume the general parameterization 



1 J. ^'_'\SHp-i'raR(t) 

7e(fc,/)(5,t) « ' + ^ ) ^ - ccnif)) (a')'"'" (a's)""^*^ 

as suggested in Ref. [28] . The trajectory is given by 

OLR^i) = 8r-\-ol {t- m%) = aR{0) + a' t , 



(12) 



(13) 



(14) 



with a' = 0.91 GeV ^. We notice the Regge poles at In — «/?(*) = 0,— 1,— 2,---. The parameters we use are reported 



Trajectory TZ 


SR 




aij(O) 


P 


1 


1 


0.5 


a2 


2 


1 


0.5 










« 



TABLE IV: Parameters of the Regge trajectories. Exchange degeneracy is assumed, 
in Table IV. Near t = m\, Eq. (13) reduces to 



{t - ml) 



(15) 



We write = /3f^/3^ using factorization of the residues at the two (1 and 2) vertices. Therefore Eq. (15) allows 
to identify as the product of two on-shell coupling constants. The residues can be obtained by the decay rates 
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p — > TTTT, ai.2 — > pTT. More precisely we obtain /3^+^o = 8.2 and /3^j^ ~ 2 from the p tttt and 04 pi: decay widths, 
respectively. Due to small value of the residue /?^^^ we will neglect the contribution of this channel in the sequel. 

Let us now discuss the 02 exchange. The residue /J^+^o can be derived from the strong coupling constant defined 
by 

M{at{p, V) - e)n°{q)) = -|^<-g^e.a;3Ae*V9^ • (16) 

'''■02 

From a2 — > pn we get ga ~ 25 GeV~^. To compute the residue we note that the a2-exchange can only occur when 
the p intermediate particles have transverse polarization. Its residue is related to ga by /J^+^o = ga/Va' ~ 13.1. 
This phenomenological value is smaller than the theoretical value given in [28] , on the basis of the Gell-Mann, Sharp 
and Wagner model [29] for the u) ^ 3n decay. In view of the theoretical uncertainties arising from the hypothesis 
of exchange-degeneracy and from the procedure we have described, we will let this parameter vary with a spread of 
±50% around a central value, i.e. we assume 

/3;^^o = 13.1 x(l± 0.50). (17) 

We note that, though the bare amplitudes B ^ pp with transversely polarized p's are suppressed (see Table I), they 
participate nevertheless in the rescattering process due to the large residue of the 02 trajectory to the p and tt. 

As shown by table IV the tt trajectory is exponentally suppressed due to « 0. Similarly we note that we have 
also computed the parameters of the ai Regge pole, but we omit this trajectory from the analysis because its intercept 
is large and negative {a{Qi) « —0.37). 

B. Charming penguins 

Charming penguins are diagrams describing the rescattering of two charmed mesons to produce two light mesons. 
Treating them as Fcynman diagrams produces a huge theoretical uncertainty. In fact to compute them one should 
employ the chiral effective theory for light and heavy mesons. However this approach cannot be extended to hard 
meson momenta and one is forced to introduce a cut-off [11-14]. To avoid the arbitrariness of this procedure one can 
describe this class of FSI by charmed Regge trajectories. This approach was followed in [19] for B pp, K*p, K*(j) 
decays and can be easily extended to B ^ tttt; we refer to this paper for details. Let us only write down the expression 
of trajectories aD{t) and aD*{t). We use Eq. (14) with sd = 0, sd* = 1 and [19] 

ao = -1.8, a' = (0.39 ± 0.12) GeV"^ , (18) 

which shows that the intercept of these trajectories is negative. Also the residues can be computed following the 
procedure of previous subsection, i.e. using the strong coupling constants go'DTv and go'D'-n (for the values of these 

constants we follow [11]). 

Differently from the case of Kt: or the K* p final states, the bare B £)(*)£)(*) amplitudes have no CKM enliance- 
mcnt. Therefore the situation is similar to the study of the B ^ pp channel in [19] where we found that, for the pp 
final state, charming penguins are less relevant than, for example, B K* p. We have checked numerically that the 
negative intercept produces a negligible contribution from the D^*^D^*^ intermediate states to S — > tttt. 

IV. NUMERICAL RESULTS AND DISCUSSION 

We present our results by taking the 7 angle as a parameter and allowing /3^+^o to vary in the range in Eq. (17). 
We compute in Fig. 1 the branching ratios B(B° — > 7r°7r°), B{B° — > tt+tt") and B{B+ — > 7r+7r°). A survey of the 
experimental results is in table V. Here we have also reported our results for Pp+^o in the range of values given in 
Eq. (17). We sec that the role of FSI is especially important for the B tt^tt" channel. 

Since we do not pretend to have presented a complete discussion of final state interactions, our result should be 
interpreted as an indication of the relevant role played by the rescattering effects when the bare amplitudes are for 
some reason small. This is confirmed by the results for the other two channels, where charged current hamiltonian is 
involved and therefore FSI play a less relevant role. Nevertheless also for the B^ — > tt+tt" channel we can see that 
FSI contribution produce a better agreement with the data. 
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Process 


B (without FSI) 


B (with FSI) 


B (cxp.) 




0.08 


0.10 - 0.65 


1.17 ±0.32 ±0.10 [30, 31] 
[32, 33] 
1.45 ± 0.29 [20] 


B° ^ TT+TT- 


8.1 


3.8 - 4.4 


4.7 ±0.6 ±0.2 [30, 31] 
4.4 ±0.6 ±0.3 [32, 33] 
4.5 ± 0.4 [20] 


B+ ^ 7r+7r° 


5.0 


3.6 - 5.0 


5.8 ±0.6 ±0.4 [30, 31] 
5.0 ±1.2 ±0.5 [32, 33] 
5.5 ± 0.6 [20] 



TABLE V: Theoretical branching ratios for B —> tttt decay channels with and without final state interactions and their 
comparison with experimental data. The column FSI is computed with P^+^o in the range given in Eq. (17). Units 10~®. 



Ir 
0.8 
0.6 
0.4 
0.2 



30 60 90 120 150 16 



5 

4.5- 
4 ■ 

3.5 
3 

2.5 



30 60 90 120 150 180 



6r 

5.5 

5 

4.5- 
4 ■ 
3.5 

3 



30 60 90 120 150 180 



30 60 90 120 150 180 30 60 90 120 150 180 30 60 90 120 150 180 

y (degrees) y (degrees) y (degrees) 



FIG. 1: Branching ratios (units 10~®) for B — » tttt as functions of the angle 7 (degrees). From left to right the decays 
-B° — > tt^tt", TT+TT^ and B^ tt+tt". Dashed, dotted, and dot-dashed lines refer to the upper, central and lower values 

of the Regge residue in Eq. (17). 



We have also computed the integrated asymmetries 

T{B° ^ 7r°7r°) - r{B" ^ 7r°7r°) 
r(BO ^ 7v%°) + r(BO ^ ttOttO) ' 
r(ffl^^±^-)-r(B°^7r+7r-) 

r(B0^7r+7r-)±r(B0^7r+7r-) ' ^ ' 

T{B- TT-^Q) - T{B+ ^ TT+TT°) 

T{B- tt-ttO) ± T{B+ tt+ttO) ' 

The results are reported in Fig. 2. For ^oo the HFAG group reports the average of the BaBar and Belle Collaborations 
as follows [20]: ^oo = 0.28 ± 0.39 . For 7 ~ 60° our result is compatible, within error with the experiment. 

Let us finally compare our results with other approaches. The exclusive B — > tttt transitions can be studied, starting 
from first principles, in the Q CD factorization approach (BBNS) [2, 3]. In this framework, all the charmless two body 
decays of B mesons have amplitudes which are shown to factorize at lowest order in 1/mf,. In other words, neglecting 
terms suppressed by heavy quark mass, the QCD factorization predicts the naive factorization ansatz [34]. In this 
framework the branching ratios for the B ^ ww decay modes are rather sensitive to the computational scheme of 
the relevant form factor. In particular, as discussed in [2], results for the tt^tt" final state depend on a parameter At 
whose precise value is unknown, but a branching ratio of the order of 10~^ could be reached. In the phenomenological 
studies [35, 36] of these processes, the authors take into account the power-suppressed and partially unknown weak 
annihilation contributions. In particular, the last fit to charmless strangclcss final state alone in the BBNS approach, 
performed in Ref. [36], reproduces the available experimental data. Our method complements the BBNS approach 
as it takes into account in a systematic way part of the power suppressed contributions, i.e. those arising from the 
final state interactions. 

Agreement with experimental data on. B ^ tttt is also obtained in [8], where the Soft CoUinear Effective Theory 



■A.00 = 
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A-Q = 
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FIG. 2: Time integrated asymmetries as defined in Eqns. (19) as functions of the angle 7 (degrees). From the left to right .4oo, 
A-{ — and ^-0. Dashed, dotted, and dot-dashed lines refer to the upper, central and lower values of the Regge residue in Eq. 
(17). 



(SCET) [4- 7] is employed. The BBNS and the SCET approaches substantially differ in treating perturbative and non- 
perturbative effects, in particular SCET predicts non negligible long-distance charming penguin contributions. We 
do not discuss the differences between QCD-factorization and the SCET as this goes beyond the limits of the present 
work. We stress however that we do not find an important role of the long distance charming penguin diagrams, but 
we find another source of long distance effects due to rescattering of the tttt and pp channels. Charming penguins 
play a minor role here because, as discussed in section III, in the Regge theory they are strongly suppressed by the 
negative intercept of the corresponding Regge trajectory. Our results are confirmed by a different analysis [12] of 
the charming penguin contributions in B — > tttt, based on an effective lagrangian approach; also in that paper these 
contributions play a lesser role the reason being, there as in the present work, the absence of any CKM enhancement 
in the bare amplitudes. 
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